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Abstract 



We are concerned with the long-time behavior of the growth-fragmen- 
tation equation. We prove fine estimates on the principal eigenfunctions 
Cu I of the growth-fragmentation operator, giving their first-order behavior 

close to and -f oo. Using these estimates we prove a spectral gap result 
by following the technique in [1], which implies that solutions decay to 
the equilibrium exponentially fast. The growth and fragmentation coef- 
ficients we consider are quite general, essentially only assumed to behave 
asymptotically like power laws. 
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1 Introduction 

In this paper we are interested in the long-time behavior of the growth-frag- 
mentation equation, commonly used as a model for cell growth and division 
and other phenomena involving fragmentation [10, 6, 4]. There are a number of 
works which study the existence and other properties of the first eigenfunctions 
(also called profiles) of the growth-fragmentation operator and its dual [7, 2, 3] 
and the convergence of solutions to equilibrium [5, 11, 9, 8, 1, 12]. These eigen- 
functions are fundamental since they give the asymptotic shape of solutions 
(i.e., the stationary solution of the rescalcd equation) and a conserved quantity 
of the time evolution. However, precise estimates on their behavior close to 
and -|-oo arc usually not given, are very rough, or arc restricted to a particular 
kind of growth or fragmentation coefficients. Our first objective is to give accu- 
rate estimates on the first eigenfunctions, valid for a wide range of growth and 
fragmentation coefficients which include most cases in which they behave like 
power laws. We give, in most cases, the first-order behavior of both first eigen- 
functions (of the growth- fragmentation operator and its dual) ; detailed results 
are given later in this introduction. 

Our second objective is to use these estimates to show that the growth- 
fragmentation operator has a spectral gap (in a certain natural Hilbert space) 
for a wide choice of the coefficients, which is interesting because it readily im- 
plies exponential convergence to equilibrium of solutions. For this we follow 
the techniques in [1], which require careful estimates on the profiles which were 
previously available only for particular growth rates (constant and linear) . Our 
results on exponential convergence to equilibrium are valid for general coef- 
ficients behaving like power laws, improving or complementing known results 
applicable to constant or linear total fragmentation rates [5, 11, 1]. However, 
our results still impose some restrictions on the fragment distribution (which 
must be bounded below) and the decay of the total fragmentation rate for small 



Let us introduce the equation under study more precisely and state our main 
results. The growth-fragmentation equation is given by 

dtgt{x) + da,{T{x)gt{x)) + \gt{x) = C[gt]{x), (la) 

(TgO(0)=0 (i>0), (lb) 

50 (a;) = .gin (a;) (x > 0). (Ic) 

The unknown is a function gt{x) which depends on the time t > and on 
a; > 0, and for which an initial condition g\-a is given at time f = 0. The 
positive function r represents the growth rate. The symbol C stands for the 
fragmentation operator (see below), and A is the largest eigenvalue of the op- 
erator g H> —dxirg) + Cg, acting on a function g = g{x) depending only on x. 
The main motivation for the study of equation (1) is the closely related 

dtnt{x) + dx{T{x)nt{x)) = C[nt]{x), (2) 

with the same initial and boundary conditions. Solutions of the two are related 
by nt{x) = e^*gt{x), and rit represents the size distribution at a given time t of 
a population of cells (or other objects) undergoing growth and division phenom- 
ena. The population grows at an exponential rate determined by A > 0, called 
the Malthus parameter, and approaches an asymptotic shape for large times. 
Equation (1) has a stationary solution and is more convenient for studying its 
asymptotic behavior, which is why it is commonly considered. Of course, results 
about (1) are easily translated to results about (2) through the simple change 
nt{x) = e^*gt{x). 

The fragmentation operator C acts on a function g = g{x) as 

Cg{x) := C+g{x) - B{x)g{x), 

where the positive part £+ is given by 

^+9{x) ■■= / b{y,x)g(y)dy. 

J X 

The coefficient b{y,x), defined for y > x > 0, is the fragmentation coefficient, 
and B{x) is the total fragmentation rate of cells of size x > 0. It is obtained 
from b through 

pX 

B{x):^ / ^b{x,y)dy (x > 0). 
Jo X 

The eigenproblem associated to (1) is the problem of finding both a station- 
ary solution and a stationary solution of the dual equation, this is, the first 
eigenfunction of the growth-fragmentation operator g i-^ ~{Tg)' + C{g) and 
of its dual If t-^ Tf' + C*{(p). If A is the largest eigenvalue of the operator 
g H> — (t g)' + Cg, the associated eigenvector G satisfies 

{t{x) G{x))' 4- \G{x) = C{G){x), (3a) 

r(.T)G(.T)|^^P = 0, (3b) 

/•OO 

G>0, / G{x)dx = l. (3c) 



Of course, the eigenvector G is an equilibrium (i.e., a stationary solution) of 
equation (1). The associated dual eigenprobleni reads 

-t{x)^'{x) + {B{x) + A) (t){x) = /:;0(x), (4a) 

0>O, / G{x)(l){x) dx = I, (4b) 

where 

C\(j){x) ;= / b{x,y)(j){y)dy. 
Jo 

This dual eigenprobleni is interesting because 4> gives a conservation law for (1): 



(j){x) gt{x) dx — I (t>{x) gin{x) dx — Cst (i > 0). 
Jo 

In this paper we always denote by G, (j> and A the solution to (3) and (4). 

In the rest of this introduction we describe the assumptions used through- 
out the paper and state our main results. In Section 2 we give the proof of 
our estimates on the stationary solution G, and Section 3 is devoted to es- 
timates of the dual eigenfunction (j>. Our results on the spectral gap of the 
growth-fragmentation operator are proved in Section 4, and we also include two 
appendices: one. Appendix A, on different approximation procedures that may 
be used for G and 0, and which are more convenient in some of our arguments; 
and Appendix B, which gives asymptotic estimates of some of the expressions 
involving the positive part £+ of the fragmentation operator, and arc used for 
our large- a; estimates of G. 

1.1 Assumptions on the coefficients 

For proving our results we need some or all of the following assumptions. First of 
all, we assume that the fragmentation coefficient b is of self-similar form, which 
is general enough to encompass most interesting examples and still allows us to 
obtain accurate results on the asymptotics of G and 0: 

Hypothesis 1.1 (Self-similar fragmentation rate). The coefficient b{x , y) is of 
the form 

b{x,y)^B{x)-p(y) (5) 

X \xy 

for some locally integrahle B : (0, -l-oo) — > (0, -|-oo), and some nonnegative finite 

measure p on [0, 1] satisfying the mass preserving condition 

1 
zp{z) dz ~ 1, 



/o 
and also the condition 

/ p{z)dz > 1. 

(When writing the integral of a measure it is always understood that the inte- 
gration limits are included in the integral.) 



The measure p gives the distribution of fragments obtained when a particle 
of a certain size breaks. 

Remark 1.1. Define, for fc > 0, the moment 

TTfe := / z'^p{z) dz. 
Jo 

We have from Hypothesis 1.1 that tti = 1 and ttq > 1. PhysicaUy, ttq represents 
the mean quantity of fragments produced by the fragmentation of one particle. 
Because of the strict inequality ttq > tti , one can deduce that p is not concen- 
trated a,t z ~ 1 {i.e. p ^ ttqSi). As a consequence we have that tt^ < 1 if fc > 1 
and TTfe > 1 if fc < 1. 

Hypothesis 1.2. The growth rate is a continuous and strictly positive function 

T : (0, +oo) -^ (0,+cx)). 

Our next assumption says that the growth rate and total fragmentation rate 
have a power-law behavior for large and small sizes: 

Hypothesis 1.3 (Asymptotics of fragmentation and drift rates). Assume that 
for some constants Bq, Boo, tq, Tqo > and 70, 7, ag, a G R 

(6) 
(7) 
(8) 
(9) 

(10) 
(11) 

to ensure the existence of a solution to the eigenproblem (see [2, 7]). 

Likewise, we impose that the distribution of small fragments behave like a 
power law: 

Hypothesis 1.4 (Behavior of p close to 0). There exist po > and /i > such 
that 

p{z) = Po z^'^^ + oiz^"-^) asz^O, (12) 

with the condition 

/i - ao + 1 > 0. (13) 

Remark 1.2. When po > condition (12) is the same as 

p{z)^Poz^-' (14) 

as 2 — >■ 0. We prefer to write it as given in order to allow for po = 0, which is 
usually not allowed in the notation (14). For instance, ii p{z) is equal to in 
a neighborhood of (such as for the mitosis case, see below), then (12) holds 
with Po = 0, but (14) does not make sense. 



B{x) 


- Bo x^" as a; -> 0, 


B{x)^ 


■^ Boo x'^ as X -^ +00, 


t{x) 


^ To a;"" as a; — > 0, 


t{x) - 


■^ Too x" as a; — > +00. 


also impose the conditions 






70 - ao + 1 > 0, 




7-a+ 1 > 0, 



To find the asymptotic behavior of the function G when x — > oo the following 
hypothesis will also be needed. 

Hypothesis 1.5 (Asymptotics to second order). Assume that t is a C^ function 
and that, for some 5 > and v > ~1, 

B{x) = Boo x'^ + 0(1-''"*') asx^ +oo, (15) 

t{x) = Too Cc" + Oix"^'^) aS X -^ +CX), (16) 

p{z)^Pi{l-zY + 0{{z-lY+^) asz^l. (17) 

Finally, to prove the entropy-entropy dissipation inequality, we will need an 
additional restriction on the fragmentation coefficient. It essentially says that p 
is uniformly bounded below by some constant p > 0, and that it behaves like a 
constant at the endpoints and 1: 
Hypothesis 1.6. There exist positive constants p, po, pi > such that 

Vz G (0, 1), p{z) > p (in the sense of measures), 

P{z) ~> _po, P{z) -¥ pi, 

and ao < 2 (which is nothing but condition (13) in the case fj. = 1). 

The reader may check that [2. Theorem 1], which gives existence and unique- 
ness of G, (j), A satisfying (3) and (4) is applicable under Hypotheses 1.1-1.4. 
We assume at least these hypotheses throughout the paper in order to ensure 
the existence of profiles. 

Let us give some common examples of coefficients satisfying the above as- 
sumptions: 

Pov^rer coefficients If we set 

b{x, y) = 2x'^^^ for x > y > 0, t{x) = x°' for x > 0, 

then all our hypotheses are satisfied when 7 — a+1 > and a < 2. Observe 
that in this case B{x) = x"^ and p{z) = 1, which satisfies Hypotheses 1.1, 
1.4 with n = I and i' = 0, and also 1.6. Since r(x) is a power, it satisfies 
Hypothesis 1.2. Hypotheses 1.3 and 1.5 are also satisfied. 

Self-similar fragmentation The previous case with r(x) = x is referred to 
as the self-similar fragmentation equation. It is closely related to the 
fragmentation equation dtgt = ^(gt) (see [3, 1]). 

Mitosis Cellular division by equal mitosis is modeled by a distribution of frag- 
ments p concentrated at a size equal to one half: 

p{z) = 2S,^.. 

This measure p satisfies Hypothesis 1.4 with po ^ pi ~ (the value of fi, 
V being irrelevant). In order to make the theory work, one has to choose 
B and r such that the rest of Hypotheses are satisfied. For instance, 
B{x) ~ x^ and t{x) = x°' with 7 — a + 1 > (and then defining b(x,y) 
through (5)) are valid choices for the same reasons as before. 



1.2 Summary of main results 

Estimates on the profiles. We describe the asymptotics of the profile G 
and give accurate bounds on the eigenvector 0. Define 



and 



^■■={ 



Px if 7 > and v = Q, 

Pij^ if7 = ,. = 0, 

if 7 > and i^ > 0, or 7 < and v > —1 



7+1 — a 
1-a ' 



where the parameters are the ones appearing in the previous hypotheses. In 
Section 2.2 we prove the following result, which improves previous estimates of 
the profile G given in [3, 1, 2] 

Theorem 1.7. Assume Hypotheses 1.1-1.5. There exists C > such that 

G{x) - Ce-'^^^'a;^-". (18) 

x—¥+oo 

This result works for all the examples given before. For all of them, it 
shows that the profile G decays exponentially for large sizes, with a precise 
exponential rate given by A(x). We observe that A(a;) behaves like 2;'^+""+^ 
(with 7+ := max{7, 0}), which is always a positive power of x. There are some 
observations about this that match intuition: the equilibrium profile decays 
faster when the total fragmentation rate is stronger for large sizes, and it decays 
slower when the growth rate is larger for large sizes. Also, it is interesting to 
notice that A does not depend on the fragment distribution (this is, p), but only 
on the total fragmentation rate B. 

The additional power x^~" which gives a correction to the exponential be- 
havior, in turn, depends only on the behavior of the distribution of fragments 
p{z) close to z = 1, this is, on fragments of size close to the size of the par- 
ticle that breaks. In the mitosis case, for example, ^ = since we obtain no 
fragments of similar size when a particle breaks. 

The behavior of G{x) for x close to depends on the power ao from Hy- 
pothesis 1.3 and the distribution of small fragments that result when a particle 
breaks. The following result is proven in Section 2.3: 

Theorem 1.8. Assume Hypotheses 1.1-1. 4 with pq > 0. If ao < 1, there exists 
C > such that 

G(x) - Gx^'-°'°. 

-(f Q^o ^ Ij there exists C > such that 

G{x) -Cx>^-\ 



This shows that G is (roughly) more concentrated close to the weaker 
the growth is for smaller sizes; and is less concentrated when there are fewer 
smaller fragments resulting from breakage. This result includes cases in which 
G{x) blows up as a; —> 0, cases in which it behaves like a constant, and cases 
in which it tends to like a power. We recall that the boundary condition is 
t{x)G{x) — ?■ as x — )■ 0, which is always ensured by /x > from Hypothesis 1.4. 

For the profile we derive the following estimates, proved in section 3, by 
the use of a maximum principle (Lemma 3.2): 

Theorem 1.9. Assume Hypotheses 1.1-1. 4- // 7 > 0, there are two positive 
constants Ci and C2 such that 

Cix< (j}{x) < C2 X, Vx > 1. (19) 

// 7 < and under the additional assumption that /i = 1 and pq > in Hypoth- 
esis 1.4, there exist two positive constants Ci and C2 such that 

Ci x^-^ < (f>{x) < C2 x''-\ Vx > 1. (20) 

Estimates of (f> are significantly harder than those of G, and they have to 
be obtained through comparison arguments. To our knowledge, this is the first 
result in which can be bounded above and below by the same power (except 
for the cases in which (f) can be found explicitly). This improves the results in 
[1] also in that it is valid for a general power-law behavior of r. 

We do not include the case 7 = in the above theorem (this is, B{x) 
asymptotic to a constant as x — > +00), but we remark that in the case of B{x) 
equal to a constant (and with the very mild condition that / 6(x, y) dy is equal to 
a constant independent of x), then (/> = 1. The case t(x) = tqx is also explicit: 
in that case, A = tq and (/)(x) = Cx for some number C > 0. 

As for the behavior at zero, we prove the following result: 

Theorem 1.10. Assume Hypotheses 1.1-1.4- Then there exists a constant 
C > such that 

</>(x) ~ Ce^(-). 

We remark that the behavior of A(x) for small x is determined by whether 
{B{x)+X)/t{x) is integrable close to x = 0. Since B{x)/t{x) is always integrable 
close to X = by hypothesis (as 70 — ao > — 1), we deduce that: 

1- If 7o < 0, then 0(x) tends to a positive constant as x — )- 0. 

2. If 7o > 0, then there are three possible cases: 

(a) If ao < I7 then again ^(x) tends to a positive constant as x — )■ 0. 

(b) If Qfo = 1, then 0(x) behaves like a positive power of x as x — )■ 0. 

(c) If ao > 1, then 0(x) decays exponentially fast as x ^> 0. 



Spectral gap. The estimates of the previous theorems allow us to prove a 
spectral gap inequality. The general relative entropy principle [8, 9] applies 
here and we have 



oo />oo 



^y cf,{x)G{x)H{u{x))dx=J J cf,{y)b{x,y)G{x) 

X iH{u{x)) - H{u{y)) + H'{u{x)){u{y) ~ u{x))) dxdy, 
where H is any function and wc denote 

In the particular case of H{x) := {x — 1)^ we define 

H[g\G] := / <pG{u-lfdx (21) 

Jo 

D[g\G]:^ / cj,{x)G{y)b{y,x){u{x) - uiy))Uydx, (22) 

Jo Jx 

and obtain that 

j^H[g\G] - -i?[g|G] < 0. 

The next result shows that H is in fact bounded by a constant times D: 

Theorem 1.11. Assume that the coejficients satisfy Hypotheses 1.1-1.6 with 
one of the following additional conditions on the exponents 70 and Uq : 

• either ao > 1, 

• or ao = 1 and 70 < 1 + A/tq, 

• or ao < 1 and 70 < 2 — ao- 

Consider also that we are in the case 7 7^ 0. Then the following inequality holds 

H[g\G] < GD[g\G], (23) 

for some constant G > and for any nonnegative measurable function g : 
(OjCxd) — > IR such that J (fig ^1. Consequently, if gt is a solution of problem (1) 
the speed of convergence to equilibrium is exponential in the L^ -weighted norm 

II • II = II • \\L^G-^,pdx), i.e., 

H[gt\G]<H[go\G]e~^' fort>0. 

Remark that in general we do not know the value of the eigenvalue A which 
appears in the assumption on 70 for the case ao = 1- Nevertheless in the case 
of the self-similar fragmentation equation {i.e. t{x) = tqx) we know by inte- 
gration of equation (3a) multiplied by x that A = tq and the condition on 70 



becomes 70 < 2. Thus Theorem 1.11 includes the result of the first part of [1, 
Theorem 1.9]. 

The main restrictions on the coefficients needed for Theorem 1.11 to hold are 
the following. First, we require Hypothesis 1.6, which says that the fragment 
distribution p should be bounded below. Consequently, this does not include 
the mitosis case and other cases in which the fragment distribution has "gaps" ; 
we refer to [5] for a proof that exponential decay does hold in that case, at least 
for a constant total fragmentation rate. Second, the exponent 70 cannot be 
too large in order to ensure that the term h{x, y) which appears in the entropy 
dissipation is not too small and can be bounded below by our methods. (An 
exception to this is the case ao > 1: in this case (j){x) decays exponentially fast 
as a; — >■ 0, and this allows us to remove the upper bound on 70.) This restriction 
on 7o might be a shortcoming of the arguments we are using; we do not know 
if there is a spectral gap when it is removed. 

On the other hand, it is remarkable that Theorem 1.11 does not place any 
restrictions on the behavior of the fragmentation or growth coefficients for large 
sizes. This is a significant improvement over [1], where the behavior at and 
+CX) of the coefficients was taken to be the same power of x, and results were 
restricted to the cases in which r is constant or linear. 

2 Estimates of the profile G 

2.1 Estimates of the moments of G 

When Hypothesis 1.3 is satisfied, we define 

C:= lim ,^-..^(^ = ^h^^l±^^h^^lPB^ i]ll (24) 
where 74- ~ niax{0,7}. Remark that, for 7 > 0, wc have the relation 

e = Pi— r^ (25) 

Lemma 2.1. Assume Hypotheses 1.1-1. 4-. For any m > I + £^ it holds that 

G(x) e^(") x"-™ dx < +00. 

Proof. As usual, wc carry out a priori estimates which can be rigorously justified 
by an approximation procedure (such as the truncated equation (50)). As G 
is integrable, it is enough to prove the convergence of the above integral on 
(xo, +cxd) for a sufficiently large xq > 0. Hence, take any xq > 0, multiply (3a) 



10 



by x^ rn^A{x) ^[f^\^ 77T, > 1 + ^ and integrate on (.tq, H-oo) to obtain 
- G(xo)e^("»)T(x-o)4-™ + (m - 1) / G(x)e^("V(a;).T-'" da; 



fV 

G{y) / e^^'''>x^-"'b{y,x)dxdy (26) 



where we have done an integration by parts on the last term. 

We first consider the case ^ > (this is, 7 > and ly = 0). From Equation (9) 
we have that for any e > there exists a;o > such that 



(m - 1) / G(a;)e^(^)T(x)x~™ dx 

J XQ 

/•OO 

> (m - 1)(1 - e)roo / G(x)e'^(^)a;"-™ dx, (27) 

and, applying Lemma B.3, also such that 

rv 
G{y) I e^^'='>x^-'"b{y,x)dxdy 

<{l + e)B^p,C' G{y)e''^y^y"~"'dy (28) 



Xo Jxa 



(observe that we have used 7 > and ly = here). Using (27) and (28) we 
obtain from (26) that 

((m - 1)(1 - e)Too - (1 + e)BooPir') /" G(x)e^(^)x"-'" dx 

J XQ 

< G(xo)e^(-")T(xo)xJ-". 

When (m — 1)(1 — e)Too — {l + e}BooPiC^^ > this gives a bound for the integral 
on the left hand side. If ttt, > 1 + ^ we can always choose e small enough for this 
to be true, because of relation (25), and it proves the result. 

The remaining case is £ = 0, this is, i^ > — 1 + ^''^ ~" . In this case we have 
to substitute (28) by the following, according to Lemma B.3: 

/ G{y) / e^^''^x^-'^h{y,x)dxdy 

Jxo Jxq 

< (1 + e)BooPiC"^"T(l + iy) G'(y)e^(^y'2/i-"+'^-('^+-"+i)(i+'') dy. (29) 

J Xo 

Since z^ > - 1 + 'y+'^-°' we have 

7-I- + 1 — a ' 

1 — rn + 7 — (7+ — a + 1)(1 + J^) < —m + a. 



11 



Thus the exponent of y on the right hand side of (29) is strictly smaller than 
a — m, so we can always find xq large enough so that 

G{y) e^(^)a;i"™6(2/,a;)da;dy <e / G{y)e^^y^y'^~'^ dy. 

Xo J Xo J Xo 

Using this and (27) in (26) wc may follow a similar reasoning as before to obtain 
the result. D 

2.2 Asymptotic estimates of G as x — )■ +cxd 

In this section we prove Theorem 1.7 by using the moment estimates in Section 
2.1. 

Proof of Theorem 1.7. We divide the proof in two steps: 

Step 1: proof that the limit is finite. Again, we carry out a priori esti- 
mates on the solution which can be fully justified by using the approximation 
(50). Let us first prove that x°'^^G{x)e^^^' has a finite limit C > as a; — >■ +00, 
and later we will show that C > 0. We use equation (3a) to obtain 

(x-«T(x)G(x)e^("))' = -ex-«-V(a;)G(a;)e^(") 

/•OO 

+ x-«e^(-M h{y,x)G{y)dy. 

J X 

Let us show that the right hand side of this last expression is integrable on 
(a;o, +00) for some xq > 0. Once we have this the result is proved, since then 
x°'~^G{x)e^'^^^ must have a limit as a; — >• +00. Integrating the right hand side 
we obtain: 



-^/ x-^-^T{x)G{x)e'^^'''> dx+ a;-«e^(^M b{y,x)G{y)Aydx 

J Xq J Xq Jx 

= /" G{x)( f y-^b{x,y)e^^y^dy~^x~i~^T{x)e^^'^n dx. (30) 

J Xg \J Xo / 

We just need to show that the parenthesis is of the order of e^'^^^a;"^''^^^'^ for 
some e > 0, and then Lemma 2.1 shows that the above integral is finite. 

The case ^ > 0. Let us start considering the case ^ > (this is, 7 > and 
1/ — 0). Using Lemma B.3 

y^^b{x, y)e^(y) dy = piB^C^ x^^^^^e^^") + 0(a;-''+"-i-^)e^("), (31) 

for some e > 0. From (16) wc also have 

Ca;"«~V(a;)e^(") = r^^x'^-^-^e'^^'^^ + 0{x'^-^-^-^)e^^^\ (32) 

12 



Using (31)-(32) and the relation (25), the parenthesis in (30) is, in absolute 
value, less than Cx°'~^~^~^e^^'^^ for some constant C > 0. Hence by Lemma 
2.1 the integral in (30) is finite, and we conclude that x°'~^G{x)e^^^^ has a finite 
limit as a; — > +00 when ^ > 0. 

The case ^ = 0. In this case we have from Lemma B.3 

bix,y)e^^y^ dy - piSooC"'""r(l + i,)a;7-(7+-"+i)(i+.')gA(.)^ 



Using the same reasoning as at the end of the proof of Lemma 2.1 we have that, 
when £, = 0, 

7 - (7+ - a + 1)(1 + iy) <a-l, 

which then shows that the right hand side of (30) is finite due to Lemma 2.1. 

Step 2: proof that C > 0. In order to show that C > in (18) set F{x) := 
r(a;)G(a;)e^'-^^ and obtain the following from (3a): 

/•oo 

F'(x) = e^(^M b{y,x)G{y)dy. (33) 



In particular, F is nondecreasing, and this is enough to conclude in the case 
^ = (since then T{x)G{x)e^'-^^ must converge to a positive quantity, so the 
same must be true of x°'G{x)e^'^^''). In the case ^ > we may bound 

F'{x) = p/(-) / b{y,x)—-e-'''^y^F{y)dy 
Jx r[y) 

poo 1 

>F{x)e^^^^ 5(y,^)^e-Afa)dy, 

Jx r[y) 

which implies that 

F{x) > F{xo) exp ( / S{w) dw, 

with 

f°° 1 

Siw) := e^^-") / b{y, w)—-^^'^^^ dy. 
Jw r[y) 



Due to equation (16) we have 



1 _ 1 

t{x) TooX° 



Riix), 
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with Ri{x) = 0{x ^ ^). Using this, and due to Lemma B.3, 

: nX pX 1 

S(w)dw> / e"^'™) / 6(y,w)— -e-^^^Mydw 

Jxo Jw ny) 

ny) Jxo 



PiBooC^ r ^{y--^+R,{y))dy 

Jxo T[y) 

= ^ r-dy+ r i?3 (y) dy > C log(y ) + Ci , 

Jxo y Jxo 



with R2{y) — 0{y°' ^ '^), Rsiy) ~ 0{y ^ '^). and Ci G IR some real number. 
As a consequence, 

F{x) >F{xo)x^e^\ 

which shows that \\v[\x^+oo F{x)x~^ (which we know exists) must be strictly 
positive. This finishes the proof. D 

2.3 Asymptotic estimates of G as a; — )■ 
Proof of Theorem 1.8. Define 

We know from [2] that F(x) — >■ when a; — ?■ and more precisely that F{x) < 
C x^ . The derivative of F, as noted in (33), is 



F'{x) = e^(-n b{y,x)G{jj)dy>Q 

J X 

so F is increasing. 

Case ao < 1 In this case, A(x) -> A(0) < 0. Choose e > such that p is a 
function on [0, e) (the fact that this can be done for small enough e is implicit 
in Hypothesis 1.4), and call p* = pljo^ej. Then, from Hypothesis 1.4, 

a;^">*(-) ^Poy^"^ asx-J-O, (34) 

with the above convergence being pointwise in y. We may additionally choose 
e e (0, 1) and C > such that 



P 



{z) KCzt"-^ for all z £ (0,e). (35) 
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Now we write 

x^-M h{y,x)G{y)dy 

J X 

^ ^Giy)A^) dy + .-^ r ^GiyU^) dy 

y ^y^ Jf y ^y^ 

For the first term in the r.h.s., we use that B{y) ~ Bgy'^" and G{y) < Gy'^~"" 

2/-J-0 

(see [2]) to write 

a;i-M f b(-]g(-]p{z)— <Cx'"'+^-"" f z°'°-^'-'°-^p{z)dz 

and conclude that it tends to zero when a; — ?■ since 70 + 1 — ao > 0. For the 
second term, we use (34) and (35) to obtain by dominated convergence 



ri-M 



°°B{y) 



G{y)p. (-) dy—~,pof B{y)y-^G{y) dy 



This hmit is strictly positive and finite, since G{y) < Cy^ "° and 70 — ao > ~1- 
Finally, we have deduced that there is a positive constant C > such that 

F'(x) ~ Cx^-\ 

which by integration gives 

t{x)G{x) ^Gx^' 

and so 

G(x) - G^ - Gx''-°"'. 

a;-i-0 T(X) x-^0 

Case ao > 1 In this case we necessarily have 7 > and 

A(a;)--Ca;^""°. 
As a consequence, following a similar reasoning as for the previous case, we have 

F'(x)^Cia;''-ie-^=^"°° 
and consequently 



due to THopital's rule. This finally gives G'(x) -~ Cx^"^ D 

x-^0 
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3 Estimates of the dual eigenfunction 

3.1 Asymptotic estimates of as a: —t- 

We first give the proof of Theorem 1.10, which is rather direct: 
Proof of Theorem 1.10. Define 

This function is decreasing since it satisfies 

'''K^) Jo 

Moreover it is a positive function, so to prove Theorem 1.10 we only have to 
prove that ip is bounded at x = 0. Consider, for 77 > 0, r^ as defined in the 
approximation procedure (see (49) in Appendix A). Then denote by 0^, A,, and 
tjjr] the corresponding functions. First we know from [2] that (/)^ converges locally 
uniformly to (j) when 77 —^ 0. We have, for ry > 0, that — A^(x) = J^ , ^V dy 
is bounded at a; = and this is why it is useful to consider this regularization. 
We have for any a;o > 0, 

supV',, =^,,(0) =i/v,(a;o)+ / — TT / b{y,z)(t)r,{z)dze~^^'^yUy 
[R+ Jo T^[y) Jo 

<M^o)+ / -^ / b{y,z)(P„{z)e-^-'-'Uzdy 
Jo T-{y) Jo 

r^o 2 ry 

<ipr,{xo) + snpiJr, — TT / b{yiz)dzdy 

Jo Tvyy) Jo 

, . N , T" B{y) fy fz\ dz , 

= .^,(xo)+supv.,y^ ^^y^^y-d, 

= i\-i {xq) + TTo sup V'r, / —TT dy. 

Jo ny) 



Now, because — is integrable at x = 0, we can choose xq > such that 
= p < 1 and we obtain 

(1 - p)sVLpipr,{x) < Ipriixa) > 1p{xo)- 

n^o 
So ipri is uniformly bounded when 77 — >■ and thus the limit i^^x) is bounded. D 
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3.2 A maximum principle 

For finding the bounds on the dual eigenfunction at a: — > +00 we use comparison 
arguments, valid for each truncated problem on [0, L] (see Appendix A for details 
on the truncation). Then we pass to the limit, as the bounds we obtain are 
independent of L. The function iphix) satisfies the equation 

S4>l{x)^Q (x6(0,L)), 

where S is the operator given by 

Sw{x) := —t{x)w'{x) + {B{x) + \i^)w{x) — / b{x,y)w{y)dy, 



defined for all functions w £ W^'°°{0,L) and for x G (0, L). This operator 
satisfies 

VluG W^i'°°(0,i) s.t. u;(i) =0, / Sw{x)GL{x)dx = (36) 

Jo 

where Gl is the eigenfunction of the truncated growth-fragmentation operator. 
We recall the concept of supers olution: 

Definition 3.1. We say that w G W^'°°{0,L) is a supersolution of S on the 
interval / C (0, L) when 

Sw{x) > (.T e /). 

Maximum principles were a powerful tool for proving the existence of sub 
and supersolutions for the growth- fragmentation models as in [7, 2]. For our 
case, we recall the maximum principle given in [2] . 

Lemma 3.2 (Maximum principle for S). Assume Hypotheses 1.1-1.3. There 
exists A> Q, independent of L, such that if w is a supersolution of S on {A^L), 
w >0 on [0, A] and w{L) > then w > on [A, L\. 

Proof. Wc start from the fact w is a supersolution on (A, L) 

~t{x)w'{x) + {B{x) + Xl)w{x) - / b{x, y)w{y) dy =: f{x) > 0. 

Jq 

Testing this equation against lw<o we obtain on (A, L) 

-t{x)w'_ {x) + {B{x) + Al)u;_ [x) = l,„(a:)<o / K^, y)w{y) dy + /(a;)l^(^)<o 

Jo 

> / h{x,y)w^{y)dy + f{x)l^(.:^)<Q. 
Jo 

Extend / by zero on [0, A]. Since w^{x) = on [0,A] by assumption, the latter 
inequality holds true on (0, L) and it writes 

yx(E{0,L), Sw-{x) > f{x)l^(^)<o. 
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Testing this last inequality against Gl, we obtain using (36) 

p L pL 

Jo J A 

Because / and Gl are positive on (A,L), this is possible only if lw<o = on 
{A, L) and it ends the proof. D 

3.3 Asymptotic estimates of as a: —t- +oo 

Now we prove the results concerning the asymptotic behavior of 4'{x) when 
X — >• +00, Theorem 1.9. For these results, we still assume that Hypotheses 1.1- 

1.4 are satisfied and, in the case 7 < 0, we additionally assume that fj, = 1 and 
Po > (so that p{z) > po > 0). We recall that Hypothesis 1.3 says that -6(2;) 

behaves like a 7-power of x and t{x) like an a-power of x, with 7 + 1 — a > 0. 

Proof of Theorem 1.9. The proof is done in two cases, and each case is 
proved in two steps. In the first step we give particular supersolutions and 
prove the upper bound, and in the second one we do the corresponding for 
lower bounds. 

Case 1: 7 > 0. 

Step 1: Upper bounds. We claim that for any C > 0, there exists ^ > and 
L* > such that 

v{x) --Gx + l-x^ 

is a supersolution on [A, L] for any L > L^, provided that max(0, a — 7) < k < 1. 
First we recall that 7+I — Q;>Oby assumption, so a — 7 < 1 and we can find 
k G (a — 7, 1). Then 

Sv{x) = -T{x)iC - kx''-^) + X{Cx - x'' + 1) + (TTfc - l)Bix)x'' - (tto - l)B{x) 

and the right hand side is positive for x large enough because the dominant 
teruiis CXx+{Trk — l)B{x)x'^ ^ CXx + {nk — 1)B aox'''^''- Indeed iTk > 1 because 
fc < 1 (sec Remark 1.1) and the dominant power is 7 + fc because fc > and 

7 + fc > a. 

Now we prove that there exists C > such that 

Vx>0, ^{x) <C{l + x). 

First we can choose C such that v{x) — Cx + I — x'' is bounded below by 
a positive constant. Moreover we take an approximation 0l of such that 
4)l{L) = 0. Then, choosing K > Q large enough, we have that Kv{x) > (j){x) 
on [0, A] because (j) is bounded uniformly in L on [0, A], and Kv{L) = KGL + 
K — KL^ > for L large enough. So, using the maximum principle and the 
previous lemma, we obtain that 

Va;>0, (J3{x) <Kv{x) <G{l + x). 
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Step 2: Lower bounds. For the lower bounds wc first prove that v{x) := 
X + x^ — 1 is a subsolution for max(0, 1 — 7) < /c < 1. 
The idea is to use x^ to transform x which is a supersolution into a subsolution. 

Sv{x) = -t{x){1 + kx^-^) + \{x + x^ - 1) - (TTfc - l)B{x)x^ + (tto - l)B{x) 

where Tr^ > 1 since fc < 1. Due to Assumption (7), B{x)x^ ^ Baox''^'^ and v{x) 
is a subsolution for x large because fc > and 7 + fc > 1 . 

For 7 > 0, there exists C > such that 

Vx>0, <p{x)>C{x-l) + . 

We know that is positive, so for C small enough, C{x + x^ — 1) — ^(.t) < 
on [0, A]. Moreover, taking an approximation (pL of (f) such that 4>l{L) ~ L, 
we have Cv{L) — (j){L) < for C < 1 and L large enough. Finally we use the 
lemmas on the maximum principle and the subsolution to conclude that there 
exists C > such that 

Vx>0, (l){x)>C{x + x'' -1) 

and the result follows. 



Step 1: Upper bounds. We start by proving that for any rj > (-g-^^ — ) "* , 



v{x) = {11 + x)'^ "'^ is a supersolution. We compute 

Sv{x) = (1 - j)Tix)irj + xy^^ + (A + B(.T))(r/ + x)''-^ 



and to estimate the last term in the r.h.s. we proceed similarly as in the proof 
of Theorem 1.8. We write, for e G (0, 1), 

b{x, y){v + yy-' dy^^^ H {r^ + y)''- V(-) dy 
X Jq \x/ 

+ B{x) / {r] + zx)'<^^p{z)dz. 



Then, choosing e such that (35) is satisfied (for this we use Hypothesis 1.4), we 
obtain by dominated convergence from (34) that 

Bi^) r% , ^7-l (V\, Bix)pov'' 



[rj + yy V(-j dy 



K— S.+00 X —7 

On the other hand we have 

B{x) I {T] + zxy-^p{z)dz - x'^^^Bix) I z^-^p{z)dz. 
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Since 7 < 0, we obtain 

a;— >-+oo — 'Y 

and finally 

x^+00 \ -7 / 

because t{x) ^ Tc^x" and a — 1 < 7 < 0. So v{x) is a supersolution for x large 

Now, we claim that there exist C > and e > such that 

Va:>0, (j)ix) <C{r] + x)'<-\ 

The proof of this fact follows from the maximum principle and taking the an 
approximation 0^ of (j) such that 4>l{L) = and that v{x) is a supersolution. 

Step 2: Lower hounds. For the lower bounds we define 

for < X < e, 



'"^^' ■ \ (x - t)x''~'^ for x> t. 

Then for e < ( "I ' | , w is a subsolution. Indeed we have for x > e 

Sv{x) = t{x){x^-^ + (7 ~ 2){x - e)x'<-^) + (A + B{x)){x - e)x''-^ 



and, reasoning as in Step 1, we obtain that 

Svix) - i\-B^poC,)x^-\ 

Finally, there exist C > and e > such that 

Vx>0, (j){x)>Cx''~^{x-e)+. 

Again, choosing an approximation (j)]^ of (f> such that (J}l{L) = L, the proof 
uses the maximum principle and the fact that vix) is a subsolution. D 

4 Entropy dissipation inequality 

As it was seen in [8, 9, 5, 1] the general relative entropy principle applies to 
solutions of (1). We remind that we use the entropy iJ[(7|G] defined in (21), 
with dissipation £'[gr|G] given by (22). We recall that 

j^H[g\G] = -D[g\G] < 0. 

For the proof of the entropy inequality we will use [1, Lemma 2.2] with C,{x) = 1. 
We need to check its hypotheses. 
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Lemma 4.1. Assume that Hypotheses 1.1-1.3 and 1.6 are satisfied with 7 7^ 0. 
Given M > \ there exists K > Q and R > 1 such that the profiles <f) and G 
satisfy the relations 

< G{x) <K {x> 0), (37) 

Giy)^iy)dy<KG{x) (x > M), (38) 

Rx 

(j){y) < K(j){z) (max{2i?M, Rz] < y < 2Rz). (39) 

Proof. The bound (37) on G is true because of Theorem 1.7 and Theorem 1.8 
with /i = 1. For the second bound, we have due to THopital's rule and using 
Theorem 1.7 

G{y)(f>{y) dy<K y^+^-^e-^^y^ dy - Kx^+^-^e'-^'^'^''"! 

Rx J Rx 

< ifa;^-"e-^'^) - A'G'(x). 
Finally, (39) is a consequence of Theorem 1.9. D 

Moreover, for proving the entropy-entropy dissipation inequality, we will 
need the following bounds, similar to those required in [1, Theorem 2.4]. 

Lemma 4.2. Suppose that the eoefficients satisfy Hypotheses 1.1-1.3 and 1.6 

with one of the following additional conditions on the exponents 70 and ao .' 

• either ag > 1, 

• or ao = 1 and 70 < 1 + A/tq, 

• or Qfo < 1 and 70 < 2 — ao- 

Let G and (f> be the stationary profiles for the problems (3) and (4). Then we 
can choose constants K, M > and R> \ such that the profiles (j) and G satisfy 

• // 7 > then 

G{x)(j){y) < Kb{y, x) {0 < x < y < max{2i?a;, 2i?M}), (40) 

y-i < Kb{y, x) [y >M, y>x>0). (41) 

• // 7 < then 

G{x)(j>{y)<Kb{y,x) (0 < x < y). (42) 

Proof. Step 1: 7 > 0, y < 2RM and x < y. We need to estimate G{x)4>{y) 
at the limit x < y — > 0. Assume for the moment that ciq < I. Using Theo- 
rem 1.8 {G{x) ~ Cx^~"° , notice that due to Hypothesis 1.6 one has fJ, ~ 1) and 
Theorem 1.10 to bound G{x) and (f>{y) respectively, we have 

Gix)(Piy) < Cx^-"°e^^y^ < C' y^'"'' 
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since ao < 1. Then under the condition 70 < 2 — ao and from Hypothesis 1.6, 
we get 

G{x)^{y) < Cy''"-' < Kb{y,x). (43) 

When ao = 1, we can do better since in this case we have necessarily 70 > 
and 

A(y) ~ — in(y)- 

y-i-O To 

Thus we can write 

G{x)<P{y) < Cy^^^° 

and we obtain the bound G{x)(f>(y) < Kh{y,x) from Hypothesis 1.6 as soon as 
7o < 1 - — . 

When ofQ > 1 we know that (j){y) decays exponentially as 2/ — ?► 0, and one 
easily sees that the bound (43) holds without any restriction on 79. 

Step 2: 7 > and 2RM < y < 2Rx. We need to estimate G{x)(j){y) at the 
limit 2Rx >y^f +00. Using (19) and (18) we have 

G[x)c^{y) < C(l + y)x«-"e-^(") 

<C'(l + 2;)y«-"e-^(^/2«) 
< C"y''-^ 

where G" depends on a, 7 and R. We conclude by using Hypothesis 1.6. 

Step 3: 7 > 0, y > A/ and y > x > 0. Since /i = 1 by assumption, we know 
from Theorems 1.8 and 1.7 that G{x) is bounded. When 7 > 0, we observe first 
that y~^ < Gy^~^ and we conclude that (41) holds true by using Hypothesis 1.6. 
Step 4: 7 < 0. When 7 < we have from Theorems 1.9, 1.10 and Hypothe- 
sis 1.6 that (t>{y) < Gy^'^^ < Kh{y,x) for ah < a; < y. D 

At this point, we have all the tools to prove the entropy - entropy dissipation 
inequality. 

Proof of Theorem 1.11. From [1, Lemma 2.1] one can rewrite the entropy 
as follows 

/>oo />oo 

D2[g\G]:= / <f>{x)G{x)<f>{y)G{y){u{x) ~ u{y))' dydx = H[g\G]. (44) 

Jo Jx 

If one looks at the integrand, one realizes that D and D2 have both 4){x) and 
G{y) as a common terms. So we would like to compare and check that 

G{xmy)<Kb{y,x). (45) 

We will denote by G any constant depending on G, 0, K^ M, or R, but not on 
g. We now distinguish two cases. 
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Case 7 < 0. The relation (45) is satisfied due to (42). So we can compare 
pointwise the integrands of D2[g\G] with D[g|G] and the inequahty (23) holds. 

Case 7 > 0. For proving the case 7 > we follow the same argument as in 
[1, Theorem 2.4]. We start by rewriting Z32[<7|G] as follows: 

D2[g\G] ^ D2.i[9\G] + D2MGI 

where 

D2,^ ■■= J J cl)ix)Gix)(biy)G{y){u{x) - u{y))^ dydx 

with Ai := {{x, y) e R']_ : y > x , y < RM or y < Rx} and A2 == A^. For the 
first term and thanks to (40) one has 

D2.i[g\G]< / Kh{y,x)ct){x)G{y){u{x)-u{y)fdydx 

Jo Jx 

<KD[g\G]. (46) 

For the other term, what we have is 

D2,2[9\G]<C / y-'<j>{x)G{y){u{x)~u{y)fdydx 

Jo Jmax{x,M} 

/•OO /"OO 

<CK / b{y,x)cb{x)G{y){u{x)~u{y)fdydx 

Jo Jma,K{x,M} 

<GKD[g\G], (47) 

where in the first inequality we applied [1, Lemma 2.2] with the bounds given 
in Lemma 4.2 and for the second one we used (41). The proof concludes by 
gathering (46) and (47). D 

A Approximation procedures 

To prove the estimates on the dual eigenfunction (p, we use a truncated problem. 
More precisely, we use alternatively one of the following ones, which differ only 
in their boundary condition 

-T{x)d^q^L{x) + {B{x) + XL)(t>L{^) = ^;(0l)(2^) for X G (0,i), 
(f>L{L)^0 or (I)l{L)^S>0 or (f>L{L) = 6L, ^43^ 

^L>0, [ G(x)0L(x)dx = L 

Jo 

The following lemma ensures that these truncations converge to the accurate 
limit when L -> +00. 
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Lemma A.l. There exists Lq > such that for each L > Lq the problem (48) 
has a unique solution {\l,4'l) with Al > and (jj^ S Wi^^{R+). Moreover we 
have 

Xl ^ A, 

L— S-+00 

\fA > 0, (f>L !• uniformly on [0,A). 

L— >+oo 

Proof. Wc start with the case (I>l{L) = by following the method in [2]. Define 
for ry > 

T (x) - I '' for < X < 7?, 

Then consider for e > and L > the truncated (and regularized) eigenvalue 
problem on [0, L] 

^(t,(x)Gl(.t)) + {B{x)+Xl)Gl{x) = I h{y,x)GL{y)dy, 
t„Gl(0) = eXf GL{y) dy, Gl{x) > 0, J^ GL{x)dx = 1, 

d f^ 

~T^{x)-g^^L{x) + {B{x) + Xl)c^l{x) - / b{x,y)(f>L{y) dy ^ T,^{0)e(f>L{0), 

4>l{L) = 0, 4>l{x) > 0, /ff 0L(a;)Gi(x)d.T = 1. 

(50) 

Notice that in this problem, the eigenelements (Al, Gl, (?!>l) depend on rj, 
and e and should be denoted (A^'^, G^'*^, 0^''^). We forget here the superscripts 
for the sake of clarity. 

The existence of a solution to Problem (50) is proved in the Appendix of [2] 
by using the Krein-Rutman theorem. Then we need to pass to the limit ry, e — > 0. 
The uniform estimates in [2] allow to do so, provided that A^''^ is positive for 
all 77, e. In [2] this condition is ensured for L large enough under the constraint 
that eL is a fixed constant, which means that L = L{e) tends to +00 as e — >■ 0. 
Here we want to pass to the limit e — > for a fixed positive value of L. For this 
we prove the existence of a constant Lo > such that A^''^ > for all 7;, e > 
and all L > Lq- 

Assume by contradiction that Al < 0. Then we have by integration of the 
direct eigenequation between and x < L 



r G{y). 
Jo 



> A / G(y) dy 

-t{x)G{x)~ [ B{y)G{y)dy+ f f b{y,z)G{y)dydz 

Jo Jo Jz 

-t{x)G{x) + (^0 - 1) f B{y)G{y) dy + J^ (J' b{y, z) dz ] G{y) dz. 
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We assume that b{y,x) = ^^p(^) with /^ p{h)dh = ttq > 1. Thus, for p 
bounded, there exists s G (0, 1) such that J„ p{h) dh > ttq ~ 1. For L > y > x > 
sL, we have 



b{y,z)dz = B{y) / p{h) dh 
Jo 



> B(2/) / p(/i) dh > B{y) / p(/i) dh > (tto - l)S(y), 
Jo Jo 



so for all X > sL 



> -T(a-)G(a-) + (^0 - 1) /" i?(2/)G(y) dy 

Jo 

which leads to 

B{x)G{x) > {no - 1)^ / B{y)G{y) dy > K - 1)^ / B{y)G{y) dy 
'''[^) Jo '''y.^) JsL 

and finally, by integration on \sL, L], 

(tto-I) / ^dy<l. (51) 

JsL ny) 

We have from Hypothesis 1.3 that 

3A>0, Vx>A ^^> ' 



r{x) (7ro-l)|ln(s)| 

SO, for L > — , we obtain 

" — s ' 

JsL ny) |ln(s)| J^L y 

which contradicts (51). Finally, Al > for aU L > Lq := ^. 

We have proved the existence of solution for Problem (48) in the case 
(I)l{L) = and we know that 

\l > A, 

Gl > G in L\R+), 

L— )-+oo 

VA > 0, 0L > (f> uniformly on [0, A). 

L— >+oo 

Now we use this result to treat the cases 4>l[L) = (5 > and (J)l{Ij) = SL. Since 
(5 > 0, we can prove by using the Krein-Rutman theorem the existence of a 



25 



solution to Problem (48). To prove the convergence of A^ to A, we integrate the 
equation on (p^ multiplied by G and we obtain 

A-AL=T(L)G(L)0i(L). 

We know from estimates on G that t{L)G(L)L — >■ when L — s> +00, which 
ensures the convergence of A. Because A > 0, it also ensures the existence of Lq 
such that \l > for L > Lq, which allows to prove the convergence of 0l to (j) 
locally uniformly (sec [2] for details). D 

B Laplace's method 

Laplace's method (see [13, II. 1, Theorem 1] for example) gives a way to calculate 
the asymptotic behavior of integrals which contain an exponential term with a 
large factor in the exponent. We give here a result of this kind, with conditions 
which are adapted to the situation encountered in Section 2. 

Lemma B.l. Take xq,Dq G R. Assume that g : [a;o,+oo) — > IR is a measure 
and h : [xg, +00) x [Do, +00) — > IR a m,easurahle function satisfying 

g{x) ^ go{x — xo)'^ as x ^ xq, for some go 7^ and a > —1, (52) 

h{x,D) ~h{xo,D) r^ ha{x-xo)'^ 

(00) 
as X -^ Xq and D — !• +00, for some /io,w > 0, 

[g(x)|e-^o''("^'^) dx < Go for some Go > and all D>Do. (54) 

Assume also that for all D > Do, the function x >-¥ h{x,D) (with D fixed) 
attains its unique global minimum at x = xo, in the following strong sense: 
there exists a nondecreasing strictly positive function 9 : (0, +00) — > (0, +00) 
such that 

h(x, D) — h{xo, D) > 6{x — xq) for all x > xo and all D > Do- (55) 

Then, as D ^ +00, 

/"OO 

e-Dhix.D)gi^^^ dx - goD^^e'^''^'^"'^^ / x'^e"''"^^ dx. (56) 

xo JQ 

The constants implicit in (56) depend only on the constants implicit or explicit 
in (52) -(55). 

Some remarks on the conventions used above are in order. Although g is a 
measure we denote it as a function in the expressions in which it appears. For 
example, integrals in which g appears should be considered as integrals with 
respect to the measure g. Also, in equation (52), it is understood that close to 
Xo the measure g is equal to a function, and the asymptotic approximation (52) 
holds. 
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Proof. First of all, by translating g and h we may consider always that xq = 0. 
We may also assume that h{xo,D) = for all D > Dq, as otherwise one 
obviously obtains the additional factor e~^'^^'^^'^\ 

An important part of the argument is based on the observation that if one 
excludes a small region close to 0, then the rest of the integral decreases fast as 
D -^ +oo: from (53) and (55) we deduce that for some p > 



h{x,D) > pmin{l,a;"} for al\x>0, D> Do- 
Then, for D > Do and < e < 1 we have from (57): 



(57) 



gix)e-^'''^^'"Ux 

/OO 
\g{x)\e~''°''^-'"Ux < Coe^^^-^")"^", 

due to (54). If we take e := D^^u then for all D > Dq wc have 



D^2Z 



^ g{x)e-'^^(='^^Ux 



<CQe-V^-^) 



This quantity decreases faster than any power of Z? as D -^ +oo. 

For the remaining part of the integral, since we are integrating in a region 
which is closer and closer to it is easy to see due to (52) and (53) that for all 
e > there exists D. > Q such that 



{l~e)gox''e-^^^+'^''°^" dx< / gix)e-°''^''^^Ux 

Jo 



< 



D^2i: 



(1 + e)5oa;'"e-°(^-<^)''"^" dx (58) 



for all D > Df. Through the change of variables z = xD^^^ we see that 



D^2i: 



(1 - e)goa;"e-°(i+')''"^'' dx 



(1 - e)g^ 



Jo 



^a^-{l+e)hoz'^ dz 



(1 - e)goD' 



_,a^-{l+e)hoz- ^ 



where the '■~' sign denotes asymptotics as _D — > +oo. Carrying out a similar 
calculation for the last integral in (58) and letting e — !► we deduce our result. 

D 

For the next result we recall that 7+ = max{0,7} and C, is defined by (24). 
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Lemma B.2. Assume Hypotheses 1.1-1.5. There is e > such that 



Xo 



e 

x— >+oo ^ 



Proof. We use THopital's rule to calculate the limit as a; — ^ +00 of 

r^mp^A{x) 

Differentiating both the numerator and denominator we find that this limit is 
the same as the limit as x — > +cx3 of 

mx"^-i + :,m Bi£HA -■ Td {x) ' 

Using (15) and (16), we obtain that 

t{x) 

for some 5 > Q. Observing that 7+ — a + 1 > and calling e :— min{(5, 7+ — a + 
1} > we have 

Tn{x)=0{x''-'). 

In a similar way, 

To(x) = x™+T+-" + 0(x'"+'^+-"-*), 
so from (59) we obtain that the limit is whenever 

TO + 7+ — Q>/c — e, i.e., T?T,>fc — 7+ + Q — e. 
This shows the result. D 

We now use this to prove an estimate which is needed in Section 2: 
Lemma B.3. Assume Hypotheses 1.1-1. 4 with pi > 0, and take fc G R. Then, 

e^^y^y''b{x,y)dy - piSoor'""r(i/+l)x'=+^-<^+-"+^)(i+")e^("^ (60) 

x—¥-\-oa 

If we also assume Hypothesis 1.5 and v = Q (and now we allow any pi > 0) 
then there is e > such that 

e^'^y^yHix,y)dy 

Xo 

^PlBooC^X^+^-^++°'^^e^^^'^ + O (a;fe+7-7++"-l-e)gA(:.)^ ^g^) 
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XO 



Xo 



Proof. We call p* and p* , respectively, the parts of the measure p on the intervals 
[0, 1/2) and [1/2, 1], i.e., p* := pl[o,i/2) and p* := pl[i/2,i]. With this we break 
the integral we want to estimate in two parts: 

I{x):^ l\^(y)y^p{l^ dy 

Mv)y'^p, (l^ dy + y" e^^y^y^p* (^) dy =: ^(x) + r{x). 

The first part, /*, can be estimated by 

< e^'^'=/^^max{x'',x^} f p, (-) dy < iroxe^^'^^^^ ma^{x'',x^o}- 

Jxo ^■^ 

Since we will show that I*{x) behaves as given in the statement, this term is of 
lower order (since A(x) is asymptotic to a positive power of a; as a; ^- +oo) and 
can be disregarded. 

For /* wc make the change of variables z — y/x and denote D := 2;7+-a+i 
to obtain 

I*{x) := r e^^y^y^p* (I) dy = x''+^ f e^("")zV(z) dz 

= a;'^+i / exp{-Dh{z,D))g{z)dz (62) 

•^-ax{^,i} 

with 

h{z,D) := -^A{xz) = -1a (i^^T^^z) , g{z) := z>^p{z) . 

Now, the asymptotics in D of the integral in (62) can be obtained from Lemma 
B.l with aiQ = 1. Let us see that h and g indeed satisiy the needed hypotheses. 
The property (52) is satisfied with go = Pi and a = v due to Hypothesis 1.4, and 
to show (53) we write (with asymptotics notation understood to be for z ^' 1 

and D — > +cx)) 

h{z, D) - /i(l, i?) = 1 (A(x-) - A {xz)) = 1 r h±^ du 
D D J^^ t{u) 

(63) 
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which corresponds to ho = (^, u = 1 in Lemma B.l. For (54) we write 
exp {—Doh{z, D)) g{z) dz 

- / '^^P ( Tr^(^-^) ) ■^^(■z) dz 



Do 



e 



< exp ( -jA(a;) 1 / z''p{z)dz<Co 

for some Co > (which in particular depends on fc), since x n> A{x)/D = 
A(a;)/x^+^"+^ is bomided for a; > 1. This gives (54). Obviously z i— ?> h{z,D) 
attains its minimum at z = 1, and (55) is a consequence of (63) and the fact 
that h{z,D) — h{l,D) is decreasing in z for all D. 
We may then apply Lemma B.l to obtain 

/•oo 

Jq 

^PiC'^'Til + ^)2;fe+i-(7+-"+i){i+-)eA(-). 

Since /*(x) was shown to be of lower order, this is enough to show (60). 
Finally, in order to show (61), wc have 

2^(^)/p (I) '^y = f '-^'^'^y' {p (I) - p^) '^y + p^ f ^""^'^y' '^y- 

For the first term, using (17) and (60) we have 

Mv)yk [p (^) _p,) dy = o(x'=+i-(^+-"+i)(i+^))e^(-), 
and for the second term. Lemma B.2 gives 

e^^y^y^ dy = (^-^x^-i++o^e^{^) + 0{x^-^++"~^)e^^''\ 
Since 7_|. — a + 1 > 0, this shows the result. □ 
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